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Abstract 



It is well known, although sometimes ignored, that not only the d = 5 but 
also d = 6 proton decay depends on fermion mixings. In general we study 
carefully the dependence of d = 6 decay on fermion mixings using the effective 
operator approach. We find that without specifying a theory it is impossible 
to make clear predictions. Even in a given model, it is often not possible to 
determine all the physical parameters. We point out that it is possible to make 
a clear test of any grand unified theory with symmetric Yukawa couplings. 
We discuss in some detail realistic theories based on SU{5) and SO{W) gauge 
symmetry. 

I. INTRODUCTION 

The decay of the proton is the most dramatic prediction coming from matter unification. 
Since the paper by Pati and Salam in 1973 [1], proton decay has been the most important 
constraint for grand unified theories [2-13]. There are different operators contributing to the 
nucleon decay in GUTs, in supersymmetric scenarios the d = 5 contributions are the most 
important, but quite model dependent. They depend on the whole SUSY spectrum, on the 
structure of the Higgs sector and on fermion masses. In recent years these contributions 
have been under discussion, in order to understand if the minimal supersymmetric SU{5) 
[2,3] is ruled out [14,15]. There are several solutions to this very important issue in the 
context of the minimal SUSY SU(5) [16,17]. 

The d = 6 contributions for proton decay in general are the second more important, but 
they are less model dependent. From the non-diagonal part of the gauge field we get the 
gauge contributions, which basically depend only on fermion masses. The remaining d = 6 



operators coming from the Higgs sector are less important and they are quite model depen- 
dent, since we can have different structures in the Higgs sector. There are several models 
where due to a specific structure of the ffiggs sector, the d = 5 operators contributing to 
the decay of the proton are not present [18]. 

In general we study in detail the gauge d — 6 contributions. Assuming that in the fu- 
ture the decay of the proton will be measured, we analyze all possible information that we 
could get from these experiments. Using this information we will study the possibility to 
test the realistic grand unified theories based on the SU{5) and SO (10) gauge groups. Our 
analysis is valid in supersymmetric and non-supersymmetric GUT scenarios. 



II. D=6 OPERATORS 

Using the properties of the Standard Model fields we can write down the possible d = 6 
operators contributing to the decay of the proton, which are SU{3)c x SU{2)l x U{1)y 
invariant [5-7]: 

Of~^ = kf eijk ea/3 ug Qjaa 7^ Qkfib (1) 

Off^ = kj eijk ea/3 ug 7^ Qjaa dg 7^ Lf^b (2) 
Ofn^ = kl eijk dg^ 7^ Qj^a ug^ 7^ Lab (3) 

Ofy^ = kl eijk ea/3 Wa 1" Qjpa W ln Qkab (4) 

In the above expressions ki — qgutM'^-^y)-! ^"^^ ^2 = QgutM'^, y,y where M(x,y), -^(x',y') 
~ Mqjjt ~ 10^^ GeV and Qgut are the masses of the superheavy gauge bosons and the 
coupling at the GUT scale. Q = {u,d), L — {u, e); i, j and k are the color indices, a and b 
are the family indices, and a, (3 = 1,2. 

The effective operators Of~^ and Off^ (eqs. 1 and 2) appear when we integrate out the 
superheavy gauge fields {X, Y) — (3, 2, 5/3), where the X and Y fields have electric charge 
4/3 and 1/3 respectively. This is the case in theories based on the gauge group SU(5). 
Integrating out {X',Y') = (3,2, —1/3) we obtain the operators Offj^ and Ofy^ (eqs. 3 and 
4), the electric charge of Y is —2/3, while X has the same charge as Y. In 5*0(10) theories 
all these superheavy fields are present. Notice that all these operators conserve B — L, i.e. 
the proton always decays into an antilepton. A second selection rule AS"/ AS = —1,0 is 
satisfied for those operators [19]. 

Using the operators listed above, we can write the effective operators for each decay channel 
in the physical basis: 

0{,eg, dp) = kl c{eg, dp) eijk ug Y Uj eg 7^ dkp (5) 

0(e«, dg) = c{ea, dg) eijk ug 7^" Uj dgp 7^, ta (6) 

0{vu da, d^) = c{ui, da, d^p) eijk ug 7^ dja dgg 7^ ui (7) 

0{ug, da, dp) = kl c{ug, da, dp) eijk dgp 7^" Uj ug 7^ dka (8) 
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where: 



c{eZ dp) = V,''vf + {V.V^uY^V.Vi^r^ (9) 

c(e„, <fp) = kl V^^vt + ^2' {VAVij,f\V,VuDVlv^f^ (10) 

c(i/,, (^^) = A;^ (14T^c/i^)'"(T^3T4iv)^' + A;^' Vt{VrVuDVlv^VENf' (11) 

c(j.p, d^) = {V,V^^f\ui^V,r + l^r(f/l;ivV^2VflB)"; a = /3 ^ 2 (12) 

The mixing matrices Vi = U^U, V2 = E]jD, = D^^c^, = D^D, Vud = UW, Ven = 
E^N and Uen = E^^N^. The quark mixings are given by Vud = U^D = KiVckmK2, 
where Ki and K2 are diagonal matrices containing three and two phases respectively. The 
leptonic mixing Ven = K^iVi^K^ in case of Dirac neutrino, or Ven = K^Vi^ in the Majorana 
case, Vi^ and VJ-^ are the leptonic mixings at low energy in the Dirac and Majorana case 
respectively. 

Notice that in general to predict the lifetime of the proton due to the presence of 0? = 6 op- 
erators we have to know ki, k2, V^^, V2, V3, V4 and Uen- In addition we have three diagonal 
matrices containing CP violation phases, Ki, K2 and K3, in the case that the neutrino is 
Majorana. In the Dirac case there is an extra matrix with two more phases. 



III. TWO BODIES DECAY CHANNELS OF THE NUCLEON 



As we know the gauge d = 6 operators conserve B — L, therefore the nucleon decays 
into a meson and an antilepton. Let us analyze all different channels. Assuming that in the 
proton decay experiments [20] one can not distinguish the flavour of the neutrino and the 
chirality of charged leptons in the exit channel, and using the chiral Lagrangian techniques 
(see reference [21]), the decay rate of the different channels due to the presence of the gauge 
d — 6 operators are given by: 



r(p ^ K+j?) 
r(p — >• 7i~^u) 

r(p ^ r]e^) 
r(p ^ K'ep 



i=l 



c(^„ d, s^) + [1 + ^{D + 3F)]c{u,, s, d^) 



Al\a\\l + D + FfY.Uu.^d.d'^) 



(ml - 



i=l 



Al \a\^{l + D- ?,Ff{ c{ep, d^) ^ + k^ c(e^, d) ^} 



487r/>3 



^i|ar(l + L> + F)^{c(e^,d^) +ktcie'^,d) } 



r(n ^ K^V) 



X c(i.„ c^, s^)\\ + ^{D- 3F)] - c{n. s, d^)[l + ^{D + 3F)] 
^^Al\af{l + D + Frj:\c{u,,d,d^) 



i=l 



(ml — m?^'^ 



r(n ^ 7r°I7) 

r(n ^ 7r-e+) = (1 + + F)^{|c(e;3, + |c(e^, f } 



(19) 

(20) 
(21) 



In the above equations is an average Baryon mass satisfying rriB ~ m-s ~ mx, D, F 
and a are the parameters of the chiral lagrangian, and all other notation follows [21]. Here 
all coefficients of four-fermion operators are evaluated at Mz scale. takes into account 
renormalization from Mz to 1 GeV. i/i — i/g, i/^, 1/^- and = e, /i. 



Let us analyze all different channels. When the nucleon decays into a strange meson plus 
an antineutrino the amplitudes (eqs. 13 and 18) of these channels are proportional to a 
linear combination of the coefficients c(z/i, s, d'") and 0(1/1, d, s*") . In the case of the nucleon 
decays into a light unfiavored meson plus an antineutrino, the amplitudes (eqs. 14, 19 and 
20) are proportional to c{ui, d, d'^). Looking at the channels with a charged antilepton, 
we see that the amplitudes (eqs. 15, 17 and 21) of the channels with a light meson are 
proportional to a linear combination of the coefficients ci^e^^d^) and c(e^,rf), while in the 
case that we have a strange meson they are proportional to a linear combination of c(eQ,, s*") 
and c(e^, s) (eq. 16). If the neutrinos are Dirac-like we have extra channels to the decay of 
the nucleon, where we have the decays into vf and a meson. The amplitudes in this case 
are proportional to c(z/f , rf, c?*^) , c{vf, s, d^) and c(z/f , rf, s^) respectively. Notice that from 
the radiative decays [22] we get the same information as in the case of the decays into a 
charged antilepton. 

Note that from the equations 13-21 we can get only seven relations for all coefficients of the 
gauge d=6 operators contributing to nucleon decay. Therefore, if we want to test a grand 
unified theory the number of physical quantities entering in the proton decay amplitude 
must be less than seven. This is an important result which helps us to know when it is 
possible to test a GUT scenario. However, as we will see in the next section looking only at 
the antineutrino channels we can get interesting predictions. 



IV. TESTING GUT MODELS 



Let us analyze the possibility to test the realistic grand unified models, the SU (5) and 
5*0(10) theories respectively. Let us make an analysis of the operators in each theory, and 
study the physical parameters entering in the predictions for proton decay. Here we do not 
assume any particular model for fermion masses, in order to be sure that we can test the 
grand unification idea. 
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In these models the diagonahzation of the Yukawa matrices is given by: 





(22) 




(23) 


eIYeE^ y*"^ 


(24) 



• a GUT based on SU{^) 

Let us start with the simplest grand unified theories, which are based on the gauge group 
S'f/(5). In these theories the unification of quark and leptons is realized in two irreducible 
representations, 10 and 5. The minimal Higgs sector is composed by the adjoint representa- 
tion E, and two Higgses and bu in the fundamental and anti- fundamental representations 
respectively [2,3], if we want to keep the minimal Higgs sector and write down a realistic 
SU (5) theory, we need to introduce non-renormalizable operators, Planck suppressed oper- 
ators, to get the correct quark-lcpton mass relations. A second possibility is introduce a 
Higgs in the 45// representation. 

In this case we have only the operators Of~^ (eq. 1), and O^f^ (eq. 2) contributing to the 
decay of the proton. Using the equations 9-12, and taking A;2 = the coefficients for the 
proton decay predictions are given by: 



c(e^, d,)sui,) = V^'vf + {V,VuDf^V,V^j,r' (25) 

c{e^,<fp)sm,)=klV^Wt (26) 

c{uu rfJ)5t/(5) = kl iy^^uoY'^iy^VENf- a = 13^2 (27) 

O{u^,da,df)sui5)=0. (28) 



We see from these expressions that in order to make predictions in any theory based on the 
SU{5) gauge group using proton decay, we have to know ki, V^^ and the matrices V2, and 
V3. Note that in a SU(5) theory there are not decays into a z/*", even if the neutrino is a 
Dirac-like particle (see eq. 28), it could be a possibility to distinguish a SU{5) theory of the 
rest of GUTs. 

In order to compute the decay rate into antineutrinos we must use the following relation: 

ciui, da, dpYsui^.^civu d„ ds)sui5) = kt{V*KlVSKM)'''iK*r''(V^K,VcKM)'''KrS''' (29) 

1=1 

Using this expression we can see that the antineutrino channel depends on the matrices Vi 
and Kl. Since we have only three independent equations (cqs. 13, 14 and 18) for these 
channels, it is clear that we can not test a GUT theory based on SU{5). Notice that from 
the channels with charged leptons is even more difficult to get some information, due to 
the presence of the matrices V2, V3 and the elements V^^. In the naive case without all CP 
violation sources beyond Vckm we could get the information about Vi from the nucleon 
decays into antineutrinos. 
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Let us analyze a particular case, the unrealistic minimal SU{5) model, where Yu — Y^ 
and Yd = (see reference [23]), in this case we have the following relations: 

c(e^, dpY^^l^^-^^ = {KTr^ + V^^^^Kl^{K*r'^{V^^^r'] (30) 
c(ea, d'^p^uit'^'" = ^1 (31) 
c(i.,, 4, rf^)rj'(5)'""^'" = (^:)''^i''^CKM^2"^i^; « = /3 ^ 2. (32) 



Notice that in this naive GUT model, all the channels are determined by Vckm- Unfortu- 
nately it is a prediction that we lost in the case of realistic versions of SU{5). However, if 
this modification of the theory does not change the relation Yu = Y^ , we could test a SU{5) 
theory from the nucleon decays into an antineutrino (see equation 29). 



• a GUT model with symmetric Yukawa couplings. 

There are many examples of grand unified theories with symmetric Yukawa couplings. This 
is the case of SO(IO) [4] theories with two Higgses IOh and 126ff, including the minimal 
supersymmetric SO(IO) model [24,25]. 

In reference [26] has been investigated the dependence of the d — 6 gauge contributions 
on fermion mixings. They consider two different cases, the naive minimal 5*0(10), where all 
fermion masses arise from Yukawa couplings to lO^f , and the case where we have the Higgses 
IOh and 126^^. Assuming only two generations, and neglecting the possible mixings which 
appear when the neutrino mass matrix is diagonalized, they showed approximately that the 
predictions for the decay channels p ^ vr+z/ and p — * K'^l'^ do not change in the different 
models for fermion masses. At the same time, it has been showed that the predictions for 
the decays p — > K^e^, and p —>■ ij,~^tt are quite different in these two scenarios for fermion 
masses. 



In this section wc will analyze the properties of all decays in those theories, using the 
fact that the Yukawa matrices are symmetric. We will take into account the mixings of the 
third generation and all possible CP violation effects. 



In theories with symmetric Yukawa couplings we get the following relations for the mix- 
ing matrices, Uc = UK^, Dq = DK^ and Ec = EKg, where K^, and K^. are diagonal 
matrices containing three CP violating phases. In those cases Vi = K*, V2 = K*V^^, 
V3 = K^Vde and V4 = K^. Using these relations the coefficients in equations 9-12 are given 
by: 

c(e^, d,)sym = {Kf\K:r[5^^ + Vl;'j,^K!^'{K;nvl^^r]{VnEr (33) 
c(e., 4).,^ = {Kf\K:Y'[kl5^' + A:2^(^2)""(V^c^i.M)^V^i,^ir^^](l^;,'^^) (34) 
c{i^u da. d'Dsym = {K*XKl'[kl6"'^S^^ + kl6"^6^^iK:rKf]{VcKMK2y\K*,VDEVENy' (35) 

c(^f , da, d%ym = {Kf^{KimK;r^{v^^Mrs-' + s'^^iK^riv^^MT'KuhNKvLf (36) 

with a = /3 7^ 2. 
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Notice all overall phases in the different coefficients. In order to compute the decay rate into 
an antineutrino we need the following expression: 



J2 da, df3)ly^c{iyi, d^, ds) 
1=1 



sym 



IT] 



(37) 



Using the above expression, and equation 13 we find that it is possible to determine the 
factor ki = gauT/Mf^xxy- 



k,^ 



Qi 



1/4 



l^il^ I^ckmI^ + 1^2!^ I^cIcmI^j 



1/4 



where: 



Qi 



A, 



(m^ — m\;yA\ \af' 

3mE 

3m j 



-D 



(38) 

(39) 

(40) 
(41) 



Notice that we have an expression for ki. which is independent of the unknown mixing 
matrices and the CP violating phases. In other words, we find that the amphtude of the 
decay p K^V is independent of all unknown mixings and CP violating phases, this only 
depends on the factor ki. Therefore it is a possibility to test any grand unified theory with 
symmetric Yukawa matrices through this channel. 

Once we know /ci, and using the expression 14 we can find the factor /c2, solving the following 
equation: 



A/2 ~\~ 2/^2^/]^ 



CKM 



+ ki 



CKM 



CKM 



with: 



Q2 



1 + 



/C2 ^ ki 



2 87r/,^r(p ^ TT + Z/) 

rripAllaf {1 + D + Py 

{-1 + My/' 







kilVSUl'mpAllafil + D + F)^ 
Using the condition Q2 > I, we get the following relation: 

t(p K^u) 



CKM 



(42) 
(43) 

(44) 



Tip 



> 



4 


^CKM 


fil + D + Ff 


K - mm\Ai 


2 


^CKM 




A2 


2 


T/12 
^CKM 


'1 



(45) 
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It is a clear prediction of a GUT model with symmetric Yukawa couplings. 



Using the expressions 13, 14, 18, 19, and 20 we can get the following relations: 

(46) 

(47) 
(48) 



with: 



rin - 




r{p - 


^ K+y) 


T{n - 


-^ttOz/) 


t(p- 


-> TT+Z/) 


T{n - 


rfv) 


t{p - 


->■ 7r+i>) 



ml{ml - mj^f[\Ai\ 


2 


^CKMl 


\"+\ 


\M 


2 


1/12 


f] 




2 


^ckm\ 


"+\ 


l^2| 


2 


1/12 
^CKM 





2mp 

6mpmf,(l + D + F)2 



-m2)2(l-D-3F)^ 



l + ^m_3F) (49) 



Notice that using the expressions for ki and k2 (eqs. 38 and 43), and the relation between 
the different decay rates of the neutron and the proton into an antineutrino (eqs. 45-48), we 
can conclude that it is possible to make a clear test of a grand unified theory with symmetric 
Yukawa couplings. 



As we say before, there are realistic 5*0(10) theories with symmetric Yukawa couplings. 
In a 5*0(10) theory all fermions of a family live in the IGp spinor representation [4]. In this 
case the coefficients for the gauge d=6 operators are given by the equations 9-12. 



Let us analyze the most realistic and studied SO(IO) theories, where all Yukawa couplings 
are symmetric. It is the case of theories with the IOh and/or 126h Higgses [24,25,27-34]. 
We have already studied the case of GUT models with symmetric Yukawa couplings, where 
we pointed out the possibility to make a consistent check of these theories. In order to 
predict the decay rates into charged antileptons in this case, we have to know the matrices 
K2 and Vde (see eqs. 33 and 34). In those 5*0(10) theories there is a specific expression for 
the matrix Vde'- 

^V^dKYu^'Vud - (3 tan a,o + tan a,26)KY^"'' = l/^^X:y^'"^T/^E(tan a^o - tan a^2e) (50) 

In the above expressions tanaio = Viq/v^q, and tano;i26 = 'i^i26/^i26- equation 50 we see 
explicitly the relation between the different factors entering in the proton decay predictions. 



To compute the amplitude for proton decay into charged antileptons we need the following 
expression: 

[S'' + v^iMKr{K;y'{v^KMy'] E virE{vi%r (51) 

i=l 

Therefore the amplitude of the channels with charged antileptons always depend on the 
matrices K2 and Vde- Therefore it is not possible to make a clear test of the theory through 
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those channels, they are useful to distinguish between different models for fermion masses 
with symmetric Yukawa matrices. Notice that in reference [26] has been showed that the 
predictions for the decay channel p — > l~^K^ are the same in different models for fermion 
masses, however as we can appreciate from equation 51 it is not true in the general case 
when we consider all generations and the extra CP violating phases. 

V. CONCLUSIONS 

We have studied in detail the predictions coming from the gauge d = 6 operators, the less 
model dependent contributions for proton decay. Analyzing the different decay channels, we 
find that there are only seven independent equations for the coefficients involved in the two 
bodies decay channels for proton decay. In general we could say that the number of physical 
parameters involved in those predictions must be less than seven. 

We have pointed out that it is possible to make a clear test of any grand unified theory with 
symmetric Yukawa couplings through the decay of the nucleon, since in these cases the decay 
rates of the nucleon into an antineutrino are independent of the mixings matrices and the 
new sources of CP violation beyond Vckm and V;, they depend only on the factors ki and k2. 
The relations between the decays of the proton and the neutron into an antineutrino have 
been found. Notice that it is the case of realistic grand unified theories based on the 5*0(10) 
gauge group. The predictions for the decay channels with charged leptons are not the same in 
different models for fermion masses with symmetric Yukawa couplings, therefore they could 
be useful to distinguish between different models. Our results are valid in supersymmetric 
and non-supersymmetric scenarios. 
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